Spin excitations under fields in an anisotropic bond-alternating quantum S=l chain: 

contrast with Haldane spin chains 
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Inelastic neutron scattering experiments on the S=l quasi-one-dimensional bond-alternating an- 
tiferromagnet Ni(CgD24N4)(N02)C104 have been performed under magnetic fields below and above 
a critical field He at which the energy gap closes. Normal field dependece of Zeeman splitting of 
the excited triplet modes below He has been observed, but the highest mode is unusually small 
and smears out with increasing field. This can be explained by an interaction with a low- lying two 
magnon continuum at gy = it that is present in dimerized chains but absent in uniform ones. Above 
He, we find only one excited mode, in stark contrast with three massive excitations previously ob- 
served in the structurally similar Haldane-gap material NDMAP [A. Zheludev et al., Phys. Rev. B 
68, 134438 (2003)]. 

PACS numbers: 75.40.Gb, 75.10.Jm, 75.50.Ee 



Recent experimental advances allovifed studies of a new 
exciting phenomenon, namely field-induced Bose con- 
densation of magnons in gapped quantum magnets 0. 
Particularly interesting is the case of antifcrromagnetic 
S — I spin chains. For uniform chains the ground state is 
an exotic quantum spin liquid with only short-range spin 
correlations and a gap in the excitation spectrum 
Antiferromagnetic S — 1 chains that are not uniform, 
but instead feature alternating strong and weak bonds, 
are also gcipped spin liquids except at a quantum criti- 
cal point However, for sufficiently strong alternation 
their so-called dimerized ground state is qualitatively dis- 
tinct from the Haldane state /k,5^,3aJj 8]- The differences 
are significant yet subtle, and involve the breaking of the 
so-called "hidden" symmetry. This symmetry is related 
to a non-local topological "string" order parameter |9lllo| 
that is in principle not observable experimentally. For 
example, thermodynamic properties are expected to be 
almost identical for Haldane and dimerized cases. More- 
over, field- induced Bose condensation of magnons occurs 
in both types of spin chains and is expected to lead them 
into the same (in terms of preserved symmetries of the 
wavefunction) magnetized high-field state j^. Thus, at 
a first glance, for the Haldane and the dimerized spin 
chains one can expect the spin dynamics to be similar in 
a wide range of external fields. But is it indeed true? 

The recently discovered [l^ nickel chain material 
Ni(C9H24N4)(N02)C104, abbreviated to NTENP, of- 
fers an ideal experimental approach to the above prob- 
lem. NTENP features bond-alternating anti- 
ferromagnetic chains in the gapped dimerized phase, 



as shown by magnetization and ESR measuren- 
ments 0, and neutron scattering experiments [l^ 
H^ . It is structurally similar to the Haldane ma- 
terials NENP (Ni(C2H8N2)2(N02)C104) and NDMAP 
(Ni(C5Hi4N2)2N3(PF6)) whose field behavior was previ- 
ously extensively investigated I?, 1?, ItII . NTENP also 
has comparable anisotropy and energy scales, which al- 
lows a direct comparison between these systems. In the 
present study we performed neutron scattering experi- 
ments on 98% deuterated NTENP in a wide range of ap- 
plied fields. We clearly observed Zeeman splitting of the 
excited triplet states below He, a field-induced softening 
of the spin gap, and the resulting emergence of antiferro- 
magnetic long rang e order (LRO) above [3, similar 
to what is seen in NDMAP. However, certain features 
of the spin excitation spectrum in NTENP are substan- 
tially different. Notable differences in the spin correlation 
function are present already at = 0, with even more 
striking discrepancies becoming prominent at H > He- 
First we summarize the crystal and magnetic prop- 
erties of NTENP. This compound crystallizes in the tri- 
clinic system (space group PI) [T^ with lattice constants 
a=10.747(l)A, 6=9.413(2)1, c=8,789(2)l, q:=95.52(2)°, 
/3=108.98(3)° and 7=106.83(3)° The ions are 
bridged by nitrito groups along the a axis having two dif- 
ferent bond distances of 2.142(3) and 2.432(6) A. These 
chains are weakly coupled via intervening perchlorate 
counter anions. The inversion centers are situated not 
on the Ni^+ ions but on the nitrito groups, so that no 
staggered components of the magnetic moments are ex- 
pected to be retained, thus resulting in occurence of the 
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long range order above at low temperatures like in 
NDMAP HO, 113 ■ The model spin hamiltionian is written 
as: 

Ti- — y~^[JS2i-^i ■ S2i + 5JS2i ■ S2i+i — /isSigH 

i 

+D[s!n (1) 

where J is the large exchange constant, 822-1, and 
S2i+i the 5=1 spin operators, 5 the bond-alternating 
ratio, g the g tensor of Ni^+, /xb the Bohr magneton 
and H an external magnetic field, and Z? is a single ion 
anisotropy constant. From the analysis of the magnetic 
susceptibility data, the following parameter values were 
evaluated from a comparison with numerical calcula- 
tions the large exchange constant J/Ub = 54.2 K, the 
bond alternating ratio (5=0.45, the single ion anistropy 
constant D/J = 0.25 (D/fcB = 13.6 K), and g||=2.14 for 
the chain direction. A similar bond-alternating ratio was 
also evaluated in Inelastic Neutron Scattering (INS) mea- 
surements 0, Q| . Such a bond alternating ratio corre- 
sponds to the ground state in the dimer phase. This con- 
clusion was directly confirmed by magnetization, ESR 
and INS [l^l experiments on NTENP. 

Three large single crystal samples NTENP with a typ- 
ical size of 16x12x4 mm^ were prepared for the current 
study by the method described in Ref. The mosaic 
spread of the sample was typically about 1.2°. The INS 
measurements were performed on several three-axis spec- 
trometers installed at different cold and thermal neutron 
sources (INS, IN12, IN14 and IN22 at ILL-Grenoble, V2 
at HMI-Berlin, SPINS at NIST center for neutron re- 
search and LTAS at JAERI). Specific details will be re- 
ported elsewhere. The scattering plane was the (a*,c*) 
plane and the external magnetic field was applied to the 
h axis. Hence, the field direction was not perpendicu- 
lar to the chain direction (a axis) (7 = 106.83(3)°). In 
this configuration, the a* was tilted from the a axis by 
w30° and the field dependence of the (1,0,0) magnetic 
Bragg peak measured at 100 mK clearly indicated a field- 
induced ordering transition at « 11.3 T (Fig^a)). 

As a first step in our investigation we characterized the 
system at zero field. In agreement with previous studies 
|l3|. we observed three distinct branches of gap excita- 
tions that appear as resolution-limited peaks in constant- 
q scans at the ID AF zone-center. These data are shown 
in symbols in the upper panel of Fig. |5| Lines represent 
contributions of each mode, as calculated in the Single 
Mode Approximation (SMA), assuming a parabolic dis- 
persion along the chain axis and fully taking into account 
resolution effects, as was previously done in Ref. [isj . 
Each of the three observed peaks corresponds to a par- 
ticular member of the S = 1 excitation triplet that even 
at = is split by single-ion anisotropy. The three gap 
energies were estimated to be 1.06±0.01, 1.15±0.01 and 
1.96±0.01 meV, respectively. The dispersion of the lower 
modes was measured all the way to the zone-boundary, 
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FIG. 1; (a) Field dependence of the magnetic Bragg peak 
intensity at q=(l,0,0). The solid line is a fit of the experi- 
mental data to the expression a{H-Hc)^^'' with _ffc=11.35 T 
and /3c=0.316 between 11.3 T and 12 T. The f3c value is close 
to the expected 3D-XY value. The dashed line is a guide for 
the eyes, (b) Field dependence of the gap energies measured 
in NTENP. Solid lines are calculated triplet branches in a 
magnetic field and a vertical thin broken line indicates the 
critical field He- The line above is a guide for the eyes. 

as shown by a solid line in the main panel of Fig. |21 It 
is well described by the form (ftw)^ = -f cjj^^x sin^ tt/i, 
where the wave vector transfer q=(ft,,A:,/),A=l.l meV 
and LJniax — 7.2 meV. The higher-energy mode (a broken 
line) follows a similar dispersion curve with A = 2.2 meV 
and Wmax = 7.2 meV. Apart from a somewhat larger gap 
energy and bandwidth, the observed behavior appears to 
be very similar to that previously seen in NDMAP 22]. 
The first key new finding of the present study is that at 
= the relative intensity of the higher-energy member 
of the triplet is anomalously weak. This fact is consis- 
tent with the data shown in Ref. Eli , but was previously 
overlooked. Polarization effects aside, the intensity of 
all three components of the triplet in NDMAP scale as 
l/w to a very good approximation. Through application 
of the 1-st moment sum rule for 5'(q, w) [23], it can be 
shown that such behavior directly follows from the SMA, 
provided that (L'/J)^ < 1. Experimentally, for NTENP 
(where D/J is actually smaller than in NDMAP) this 
scaling is violated, and the relative intensity of the upper 
mode is lower than expected by about a factor of 3. We 
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came to this conclusion after performing measurements 
in several Brellouin zones maintaining gn = tt, to elim- 
inate the effect of polarization-dependent coefficients in 
the unpolarized neutron scattering cross section. 

The suppression of the upper mode in NTENP be- 
comes even more apparent in applied magnetic fields. 
The two lower modes in NTENP behave almost exactly 
as in NDMAP for H applied perpendicular to the spin 
chains. However, unlike in NDMAP, the upper mode 
in NTENP vanishes in relatively modest applied fields 
of if > 4 T. This is borne out in the constant-g scans 
shown in the two lower panels of Fig. |2] Thus, even 
well below He the excitation spectra of NTENP and 
NDMAP are qualitatively different. It is interesting that 
while the intensities of excitations in NTENP behave 
anomalously, the field dependencies of the correspond- 
ing gap energies at H < He are very similar to those in 
NDMAP. These field dependencies measured in NTENP 
are plotted in Fig. [Hb). As in the case of NDMAP, 
below He these data are well described by a simple- 
minded perturbation theory calculation [2^: A±{H) — 
± + {gf^^nn/^ and Ao(7J) = A„ 

where A^,, Aj, and A^ are the gaps at if = for excita- 
tions polarized along x, y and z, respectively. The best fit 
is obtained with A2=2.0±0.1 meV, A:j;=0.97±0.02 meV, 
Ay=1.13±0.02 and 5=2.09. 

The second crucial finding of this work is that at 
H > He there is only one low-energy mode at qe = tt 
in NTENP. This can be seen from Fig. 21 that shows 
constant-g scans collected below (8.5 T), around (11.5 T) 
and above(14.5 T) the critical field {He ^ 11.3 T). These 
data are in stark contrast with similar scans measured in 
NDMAP (Fig. 1 of RefH), where three distinct modes 
are visible at H > He- As discussed above, the upper 
mode in NTENP vanishes well below He- In addition, 
in NTENP the lowest-energy mode that goes soft at He 
does not reappear at higher fields, as it does in NDMAP. 
The data for NTENP at 14.5 T show only a statistically 
insignificant hint of a peak at 0.4 meV. Scans in several 
Brilloin zones confirmed that the low-energy mode is ab- 
sent at H > He in NTENP in all polarization channels. 

Comparing the behaviors of NTENP and NDMAP, we 
conclude that the spin dynamics of anisotropic Haldane 
and dimerized S=l chains are qualitatively different both 
below and above the critical field. This fact is confirmed 
by a recent Lanczos first-principles study of the dynamic 
structure factors in NTENP and NDMAP, based on ac- 
tual measured exchange and anisotropy parameters [l^ . 
The calculated dispersion relations for NTENP (solid 
symbols in Fig.|31) match our experimental data remark- 
ably well. The solid squares and circles correspond to 
peaks in the dynamic structure factors for excitations 
polarized parallel to the easy plane {S^^ and 5^^) or per- 
pendicular to it (5^^), respectively. The calculations also 
predict that at ii = the intensity of the out-of-plane 
S^-^ mode is only 20% of that for S"^^, in agreement with 
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FIG. 2: Inelastic energy scans near q= 
H=Q, 2 and 4 T. 
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FIG. 3: Measured dispersion relation of magnetic excitations 
in NTENP (open circles). Solid symbols show the positions 
of peaks in the in-plane (squares) and out-of-plane (circles) 
dynamic structure obtained by Lanczos calculations. Solid 
and broken lines are the results of fitting to the form {Hlj)^ = 
sin^ 7r/i. Inset: excitation continuum (shaded area), 
as calculated for NTENP. The lower bound (thin solid line) 
coincides with the observed position of the observed z-axis 
mode (the highest excitation mode) . The in-plane modes are 
shown by a heavy solid line. 



the neutron results. Finally, the numerical calculations 
reproduce the key features at H > He'- three modes for 
NDMAP and only one for NTENP. 

An understanding of the underlying physical mecha- 
nism emerged from a numerical size-dependence analy- 
sis (261 that can distinguish between long-lived single- 
particle excitations and peaks in multi-particle continua. 
The key lies in the violation of translational symmetry in 
a dimerized chain. This symmetry breaking makes the 
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FIG. 4: Inelastic energy scans at q=(l.l,0,1.5) and 100 mK 
for H=8.5 (< He), 11.5 (- i?c) and 14.5 T (> H^). 



wave vectors g|| = tt and gn = equivalent. As a result, 
in NTENP around q^ir there is a low-lying continuum 
that consists of paired g|| = and q|| = tt magnons. In 
the presence of sufficiently strong anisotropy the upper 
mode at q|| = tt is pushed up in energy and can en- 
ter this continuum. Once this happens, its decay into 
a pair of lower-energy excitations becomes allowed by 
energy-momentum conservation laws. The upper mode is 
then no longer a sharp (5-function of energy, but appears 
as a asymmetric finite-width peak on the lower contin- 
uum bound (Fig. |31 inset). The situation becomes pro- 
gressively worse when an external magnetic field drives 
the upper mode to still higher energies, deeper into the 
continuum. The two lower modes remain safely below 
the lower continuum bound at all times and are thus 
not affected. This scenario is reminiscent to that for 
the classical antiferromagnet TMMC jl^. In the latter, 
one magnon optic mode interacts with the two magnon 
mode to show the avoided crossing. However, the anal- 
ogy is only superficial, since, unlike in TMMC, the single- 
particle excitations in NTENP are a triplet to begin with. 

Nothing of the sort can occur in the uniform Hal- 
dane spin chains of NDMAP, where the magnon gap at 
q\\ = is much larger than that at g|| = tt |2^. The 
lowest-energy continuum states at g|| = tt are then three- 
magnon excitations, and the corresponding continuum 
gap is much larger 26] . Even in the presence of anisotropy 
all three branches are below the continuum bounds, and 
persist as a triplet of long-lived excitations [2^ . 

A low-lying continuum is also responsible for the "dis- 
appearance" of two of the three modes in NTENP at 
H > He- Numerical calculations show that, unlike be- 
low the critical field, at H > He the lower bound of the 
continuum in S^^ is actually beneath the lowest single- 



particle state. Since at any H > both the highest- and 
lowest-energy modes have mixed polarization in the (y, z) 
plane, in the high-field phase they both become subject 
to decay into multi-particle states. At H > only the 
middle mode survives in NTENP, being polarized along 
X. For this excitation all decay channels remain closed. 

In conclusion, our experimental results clearly demon- 
strate fundamental quantum-mechanical differences be- 
tween the two exotic spin-liquid phases of uniform and 
bond-alternating integral spin chains. 
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